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ABSTRACT 

We re-analyze the Chandra ACIS spectrum of the kpc-scale jet in 
PKS 0637-752 to investigate the possible low energy cutoff in the relativistic 
electron spectrum producing the non-thermal radiation in the scenario of inverse 
Compton emission off the cosmic microwave background. This was among the 
first objects targeted by the Chandra Observatory and gives a unique opportu- 
nity to study the low energy X-ray emission free of detector contamination. As 
previously reported in the literature, the spectrum can be fit by a power law, 
with the slope predicted by the radio spectrum, modified by low energy absorp- 
tion through the Galaxy as determined from the spectrum of the quasar core 
and by HI 21 cm observations. We obtain a marginally better fit with an model 
of inverse Compton emission produced by an electron population that exhibits 
a cutoff at 7min5io between about 50 and 80 (assuming F = 5). This range for 
7min is higher than has previously been assumed in broad-band spectral fits to 
the jet emission. The observed optical flux can be used to place a lower limit 
on 7min; the constraint is not very strong, but does suggest that 7inin must be 
higher than 1 to avoid overproducing the optical emission. We investigate the 
effect of uncertainties in the column density for galactic absorption as well as 
the calibration of Chandra for these early observations. Finally, we discuss the 
implication of these limits on the jet luminosity in this source. 
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Subject headings: galaxies: jets — quasars: individual (PKS 0637-752) — X-rays: 
galaxies 



INTRODUCTION 



The initial observatio ns to focus Chandra using the pres umed point source PKS 0637-752 



a quasar a t redshift 0.651 



structure (Schwartz et al. 



200G 



Chart as et al. 



Savage. Browne, fc Boltqnlll976f) . revealed a remarkable X-ray jet 



20001 ). coincident with the radio jet reported 



by lTingay et al\ (119981 ). ISchwartz et al\ (120001 ) pointed out that the X-ray emission could not 
be explained reasonably by thermal bremsstrahlung, by an extension of the radio synchrotron 
spectrum, by synchrotron self-Compton, or by inverse- Compton (IC) from any likely source 
of seed photons provided that the magnetic field was near its equipartition value. Nonethe- 
less they concluded that IC from the same electrons producing the radio emission was the 
most likely source of the X-rays in view of the similarity of the X-ray and radio surface 
brightness profile. 



Tavecchio et al\ ( 120001 ) and ICelotti. Ghisellini &: Chiabergd (120011 ) solved the dilemma 
by invoking bulk relativistic motion of the 100 kpc scale jet. Relativistic motion with bulk 
Lorentz factor F has the effect of requiring a smaller magnetic field in the rest frame for 
minimum energy conditions, and of increasing the apparent energy density of the cosmic 
microwave background (CMB) by a factor of on average as seen by the radiating elec- 
trons. Such an IC/CMB explanation subsequently has been apphed succes sfully to explain 



the X -r ay emission from pow e rful q u a.sar and PR II radio 



(2002); [Harris fc Krawczvnskil (120021): ISiemiginowska et al\ ( 



ets: e g., b v 



2002 



2003): 



Sambruna et al 



Sambruna et al. 



mm -. Marshall et al\ (120051 ): ISchwartz et al\ ( 12006[ ): ISambruna et al\ mm and ISchwartz 



mm 



Tavecchio et al\ (120001 ) pointed out that the scenario in which low energy electrons 



underwent Compton scattering on an external photon field gave an opportunity to infer 
parameters of the low energy portion of the relativistic electron distribution for the first time. 
They estimated 7min ~ 10 for the PKS 0637-752 jet, by constructing a fiducial spectral model 
which fit the bro ad-band radio, optical and X-ray flux data. Subsequent ex aminations of 
the jet spectrum ( ICelotti, Ghisellini fc Chiabergell200ll : lUchiyama et a/.ll2005l ) all converged 
on a value for '-jmin ~ 10-20. In the present work we perform a more rigorous assessment of 
7min by performing a detailed spectral fitting of the ACIS X-ray spectrum. PKS 0637-752 
is in the unique position among identified X-ray jet sources in that, because the exposures 
were done so early in the lifetime of Chandra, the low-energy X-ray spectrum is free of the 
detector contamination and associated loss of sensitivity that affect later observations. The 
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main source of uncertainty remains the low-energy calibration of ACTS, and we contrast 
results obtained from a wider bandpass (0.3-7 keV) than is usually used for imaging spectral 
analysis using Chandra with the equivalent results if the bandpass is restricted to the better- 
calibrated range of 0.55-7 keV. 

Section [2] lists the observations and data reduction for this investigation. The develop- 
ment of the custom IC/CMB spectral models is detailed in Section [3l which are then applied 
to the data in Section H] to derive limits on 7min- The interpretation of the model fits, in- 
cluding a comparison between two different prescriptions for the calculation of the IC/CMB 
spectrum and a discussion of the effect of the shape of the electron distribution at the 7min 
cut-off, is the topic of Section [51 Section [6] summarizes the paper. 



2. OBSERVATIONS AND DATA REDUCTION 

PKS 0637-752 was observed by Chandra using the ACIS-S instrument as part of the 
initial focusing period. This early in the mission, there was no loss of low-energy sensitivity 
due to the contamination layer. We thus have an opportunity with this object to push the 
X-ray spectrum to lower energies than is currently possible. 

Table [T] lists the details of the Chandra observations that were used in this investigation, 
all of which were carried out in the faint timed mode of ACIS. We selected the observations 
based on whether the quasar and jet were positioned on the back-illuminated S3 chip, which 
excludes two observations (ObsID 1059 and 1061). 

Other observations were discarded: ObsID 1057 and 62558 because the identification of 
the quasar and jet was unclear, and ObsID 1093, 1264, and 1265 since they do not appear 
to contain any useful data. Although most of the remaining observations were not at the 
best focus, the quasar and jet are clearly separated, and the jet spectrum can be extracted 
without contributions from the quasar core. 

For each observation, the evt2 file was re-extracted from the evtl file in order to apply 
the latest calibrations, as outlined in the ACIS data preparation thread on the CIAO web 
page (CALDB version 3.2.4). Since the observations were taken early in the mission, and 
the focal temperature was different from —120 °C, no time-dependent gain or CTI correction 
was used. 

Extraction regions for the jet and background were defined individually for each obser- 
vation, since in some cases the departure from optimal focus necessitated a larger extraction 
region to include all events from the jet. The extraction regions were divided into an inner 



-4- 



and outer jet, the first of which included the faint X-ray emission of the bridge between the 
quasar core and the bright radio/optical/X-ray knots, which comprised the second region. 
The regions for one of the observations (ObsID 62554) are shown in Figure [H No significant 
differences in the spectra between these two sub-regions were found, although the low num- 
ber of counts in the bridge spectrum make a detailed comparison difficult. All subsequent 
fits used only the outer jet extraction region. 

Using the CIAO tool dmextract, the counts in each observation's source and background 
extraction regions were combined into pulse invariant (PI) channels to form the individual 
spectra. A grouping algorithm was applied to have a minimum of 2 counts per bin. Weighted 
ancillary response files (WARF) and redistribution matrices (WRMF) were created using the 
tools mkwarf and mkrmf , as appropriate for extended sources, (mkacisrmf does not apply for 
these observations due to their early observation dates.) To check for background fiares, the 
light curve of the S3 chip excluding the region around the quasar and the jet was analyzed 
visually. Data during times corresponding to the few occurrences of spikes in the count rate 
were filtered out in the extraction. 

Because the results reported in this investigation depend heavily on the correct mod- 
eling of the low-energy response of acis, alternative calibration products were extracted 
to investigate the effect of uncertainties in the calibration (Grant, C. & Bautz, M., private 
communication). Sets of WARF and WRMF were constructed using — 100°C and — 120°C 
FEF files. In addition, for the — 100°C FEF calibration products, the Quantum Efficiency 
Uniformity file was varied between the fiat version NOOOl and the version N0002 that in- 
cludes information on how the quantum efficiency varies over the S3 chip. Spectral fits were 
repeated with both the nominal set of calibration files and the above alternative calibrations. 
The results on '"j^mn are unaffected to within the uncertainties introduced by the choice of IC 
scattering kernel and the shape of the 7min cutoff (see Section [3]). However, we will report 
results obtained from both a 0.3-7 keV bandpass as well as a better-calibrated 0.55-7 keV 
bandpass. 



PKS 0637-752 was also observed by the Hubble Space Telescope (ISchwartz et a/.ll2000l ). 
The area where the three optical knots are located corresponds morphologically to the ex- 
traction region used for the X-ray analysis. The total optical flux density from the knots is 
0.574/iJy at a frequency of 4.3 x 10^^ Hz. 
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3. IC/CMB MODEL 

We wish to fit the spectrum for the jet with a model of inverse Compton scattering 
of energetic electrons on the CMB. We assume that the spectrum of the entire X-ray jet is 
emitted by a spatially homogeneous population of electrons that collectively exhibit a bulk 
motion with Lorentz factor F = (1 — The random motion of the electrons (assumed 

isotropic in the jet rest frame) is described by the distribution of the Lorentz factor 7, such 
that the energy of any given electron in the jet rest frame is ^me(? {me.: electron rest mass, 
c: speed of light; primed quantities refer to the jet rest frame, unprimed to the observer 
frame; exception: 7, which always refers to the jet rest frame). 

This distribution is modeled as a power law between 7inin and 7max, with a slope s 
{n'{'y) oc 7~'^). The upper end of the distribution is chosen as a sharp cutoff; since we 
note that the high energy cutoff is not reached in the X-ray spectrum, the exact functional 
form doesn't affect the analysis as long as 7max is chosen sufficiently high. The shape of the 
low-energy cutoff, however, has a direct effect on the model spectrum. We investigate two 
options: a step function at 7min: 



^'(^\ ^ J ^0 7 * 7mm < 7 < 7max /^v 

^'^ '0 Otherwise, ^ ^ 



and a constant electron density below 7n 



r ^0 7mf„ 7<7min 
^'(7) - I K 7~' 7min < 7 < 7max (2) 
I 7inax < 7- 

Neither of them is likely to be the actual shape of the low-energy cutoff, but they 
(together with the case when 7inin is sufficiently small as to move the cutoff outside of the 
spectrum bandpass) provide three representative cases for the characterization of the electron 
distribution based on the available X-ray data. 

j'{E[, the IC/CMB emissivity in the jet rest frame at a given energy E[ and direc- 
tion can be formulated as an integral over the invident electron and photon energies and 
their directions, taking into account their respective number densities (which are themselves 
functions of the corresponding particle's energy) and the relativistic transformations of the 
energies and directions between the jet rest frame and the individual electron's rest frames. 
The scattering of CMB photons into the X-ray band is due to electrons with Lorentz factors 
small enough to make the scattering event in the electron rest frame be safely within the 
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Thomson limit. This and other assumptions allow these integrals to be evaluated in terms 
of elementary functions. We wish to compare two different approaches that differ in the way 
the photon field is treated. 



3.1. Blumenthal & Gould Approach 



The iBlumenthal fc Gouldl (Il970l ) calculation assumes an isotropic distribution of seed 
photon directions, which is strictly speaking not applicable to a jet moving relativistically 
with respect to the CMB rest frame. However, the simplicity of the resulting expression for 
the emissivity is the reason it has been used extensively in the literature. Since both the 
electron and photon fields are treated as isotropic, the emissivity does not depend on Q[ and 
can be written as 



fiE[) oc E[ 



Tcrit 



'^f{E'„E[,j)dE',dj, 



(3) 



where f{EQ,E[,'j) = 2x\nx — 2x^ + x + 1, with x = E[/ {A'j'^Eq). 7 3> 1 is assumed 
throughout. The function / is often called the kernel of IC scattering, because it describes 
the spectrum obtained from a fixed initial electron and photon energy. 7crit, which is the 
Lorentz factor of the least-energetic electrons that can contribute to the emission at E[, is 
derived as 




(4) 



A further simplification treats the incident CMB photon field as monoenergetic at an 
energy Eq = {1 + z) T a kTcMn, where a = 2.7 is the average photon energy of a thermal 
spectrum in units of kT (IFelten &: Morrisonlll966r). is the Boltzmann constant, and Tcmb = 
2.725 K is the local CMB temperature (IFixsen et a/.lll996l ). The integral over E'q in Equation 
[3] is thus eliminated. The motivation for this simplification is that features in the observed 
X-ray spectrum are expected to be much broader than the width of the thermal distribution 
of CMB photon energies. 



In the case of the electron distribution with the sharp cut off at 7n 
final expression for the emissivity is 



(Equation [T]), the 
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/"7max 

j\E[)^E[ ^-^-'f{x)d^. (5) 

Jmax{7crit,7min) 

For the other electron distribution (Equation [2]), the equivalent expression is 

f{E[)^E[ j-^-'f{x)d^, (6) 

>^7crit 

if 7crit > 7min, and 

/'7max ^7min 

fiE[) oc E[ / 7"^-' fix) d7 + E[ / 7"' f{x) d^, (7) 

>^7min '^7crit 

if 7crit < 7min- 

The observed spectral flux j{Ei) is proportional to the emissivity in the jet rest frame 
j'{E[) when the energy shift of the photons due to the motion of the jet rest frame with 
respect to the observer is taken into account. Photons emitted at an energy E[ in the 
jet rest frame are observed at an energy Ei = 6E[/{1 + z) in the observer frame, where 
6 = [r(l — is the jet Doppler factor (with 6 = cos~^ fi the observer viewing angle of 

the jet) and z is the redshift of the quasar. Since for the determination of 7min we are only 
interested in the shape of the spectrum, and not its normalization, the implementation of 
the model simply uses the jet rest frame emissivity appropriately shifted along the energy 
axis and applies a normalization factor that is not specified in detail to best fit the observed 
number of counts. 



3.2. Aharonian &z Atoyan Approach 



In contrast to iBlumenthal fc Gouldl (jl970h , the lAharonian fc AtoyanI (jl98lh approach 
treats the incident photon field as monodirectional (anti-parallel to the jet bulk velocity 
for IC/CMB) in the jet rest frame, which is an appropriate approximation for F 1. 
Because the electron distribution in the jet rest frame is assumed isotropic, the system now 
exhibits azimuthal symmetry around the propagation direction of the jet, and the angular 
dependence of the emissivity reduces to a dependence on the polar angle 9' = cos~^ fi' of the 
outgoing photon. If 7 ^ 1 is assumed as well, then, irrespective of the scattering angle in the 
electron rest frame, the direction (in the jet rest frame) of the photon after scattering is well 
approximated by the direction of the incident electron. The integration over the incident 
electron direction thus reduces to a 5-function substitution. 
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Adopting the notation of IStawarz et al\ (120051 ) . the emissivity can be written as 



(8) 



with the IC kernel, expressed using the quantities f ' = 2 (1 — /i') E'q'~^ j {nieC^) and w' 



fiE',,E[,^,fi') = l 



In this case, 7crit evaluates to 



w 



12 



2w' 



+ 



2w 



12 



2{l-w') v'{l-w') v''^{l-w'Y 



(9) 



7crit 



2mlc^ 



il-f,')E',E[ 



1/2- 



(10) 



Again, the cond i tion that 7^1 has to be satisfied. The additional condition in 
Aharonian fc AtoyanI (119811 ). that E[ ^ EL is the same simplify i ng ass umption mentioned in 
th e explanatory text to Equat ion 7.26a in lRybicki fc LightmanI (119791 ) and to Equation 2.13 



m 



Blumenthal &: Gould! (1l970l ) and is readily satisfied for IC/CMB X-rays. While the above 



expression is derived in the Klein-Nishina regime, the authors stress that the exact formula 
for the IC spectrum (without the simplifying assumptions mentioned above) is valid for any 
values of the electron and photon energies. Since the subsequent simplifications are satisfied 
in the Thomson regime also. Equation [8] is applicable to the process of X-ray generation 
via IC/CMB. In addition, f{EQ,E[,'y,fi') does not appear to require any modifications to 
mitigate potential numerical precision issues in its evaluation in the Thomson limit. 

In the case of the electron distribution with the sharp cut off at 7min (Equation [1]), the 
final expression for the emissivity is 



f{E[,f,')^E[ / ^~^''f{EiE[,^,f,')d^. (11) 

Jmax(7crit,7min) 

For the other electron distribution (Equation [2]), the equivalent expression is 

/"7max 

j\E[, /.') oc E[ / 7"^"' f{E'„ E[, 7, /i') d7, (12) 

J 7r rit 



if 7crit > 7min, and 
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j'{E[,fi') oc E[ 



7 



7 



/(E^,E;,7,/i')d7, (13) 



'Ycrit 



if 7crit < 7 



As in the iBlumenthal &: Gould (119701 ) approach, the observed spectrum j{Ei,fj,) oc 
j'{E[,fi'). fi' is related to /i and /3 through yu' = (/5 — /i) / (1 — /9yu), when the incident 
photons in the jet rest frame are travelin g in the opposite direction to the jet. Since the jet 
viewing angle is small (iLovell et a/.ll2000l ). we assume that F = 5, which implies n = jS, and 
thus /i' = 0. 



3.3. Model Implementation 



We proceed to implement custom XSPEC (jArnaudlll996l ) models that compute spectra 
based on the above formulas. The numerical integration is done using the qsimp routine 



m 



Press et ali (119921 ). For both approaches, the two adjustable model parameters are s 
(the electron distribution power law index) and 7min- The additional parameters that are 
not adjusted in the fit are 7max, F, 6, and z. 7max is kept fixed at 10^ (corresponding to a 
high-energy cut-off around 850 MeV); z = 0.651. As mentioned before, the normalization 
of the model is arbitrary, as the correspondence between the normalization and the relevant 
physical quantities like luminosity and distance is not spelled out. 



Following iTavecchio et al\ (120001 ) and ISchwarta (120071 ) . we assume F = 5 = 10, which 



sets E'q = 1.05 X 10~^ keV. From the (5-function approximation to the energy spectrum of 
IC scattering from a fixed initial electron and photon energy (where E[ = 4/3 j"^ Eq fixed), 
the Lorentz factor for electrons scattering CMB seed photons to a given observed energy 
scales as (F(5)~^/^. Under the assumption that F = 5, this reduces to a scaling by 5^^. 
The same behavior is expected for the two kernels used in the models, although this was 
not investigated in detail. Our quoted results for 7min are thus expected to depend on the 
assumed value for S in the same manner, and we report all results for 'yrnin with this scaling 
in mind. For reasonable departures from the assumed values for F and 6, an approximate 
scaling by (F5)~^/^ is expected to remain even if F = 5 is not assumed, since the dominant 
behavior of the scattering kernel in all cases includes the photon and electron energies in the 
combination E[/{'^'^ E'q) . 
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4. RESULTS 

Given that each individual spectrum has only a small number of counts over the band- 
pass of interest and was therefore grouped to have a minimum of only two counts per channel, 
is not appropriate as a fitting statistic. Instead, the C-statistic implemented in XSPEC 
was used. Unfortunately, this statistic does not allow for a goodness-of-fit test, in the way 
the reduced is commonly used. However, a visual inspection of the fits reveals a very good 
agreement between the data and the model. In any case, the important discriminant will be 
the changes in the fitting statistic with the model parameters as well as between models. All 
subsequent results were obtained from simultaneous fits to the 21 individual spectra while 
fixing the model normalization to be the same between the spectra. 



4.1. Phenomenological Fits 



We first fit the spectra in the 0.3-7 keV bandpass with a single power law, modified by 
neutral Galactic absorption. If rin is allowed to be free, a best-fit value of (9.1 ± 0.8) x 10^° 
cm~^ is returned, and the best-fitting power law energy index is a = 0.76 ± 0.04. These 
values corre spond very well with the e xpected values: the absorbing column with the value 
reported by lDickey &: LockmanI (jl990l ) (9.1 x 10^° cm~^), and the power law i ndex with both 



the r esults from earlier investigations of the C/^anG?ra spectrum of this sou rce (jChartas et al 



20001 . a = 0.85 ±0.08, ) and the radio spectral index (ISchwartz et a/.ll2000l ). The same power 
law index (to within uncertainties) is obtained by fixing nn = 9.1 x 10^° cm~^ before fitting. 
The best-fit C-statistic value is 886.3 for 958 bins. 

In the restricted 0.55-7 keV bandpass, the fit with both the power law index and the 
absorbing column free returns nn = 3.0 x 10^° cm~^ (with large uncertainties) and a = 
0.63 ± 0.05. The C-statistic at the minimum is 808.2 for 860 bins, but it is clear that the 
low value for the absorption and the resulting fit is spurious, as the restricted bandpass is 
not very sensitive to the absorption column density. If = 9.1 x 10^° cm~^is fixed before 
fitting, the power law index increases to a = 0.84 ±0.02, once more consistent with previous 
values. The minimum value of the fitting statistic is 817.2 for 860 bins. In all subsequent 
fits, the absorbing column is kept fixed at the Galactic value. 

The normalization of the power law component (before absorption) is in all cases equal 
to about (3.1±0.5) x 10"^ photons/cm^ /s/keV, corresponding to a flux density of (21±3) nJy 
at 1 keV. This is sligh tly smaller than, but probably within la of, previous measurements 
(ISchwartz et al\ 120001 ) , which could be due to differences in the extraction region or the 
background subtraction. 
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4.2. IC/CMB Model Fits 

4-2.1. Upper Limit on 7min from X-ray Spectrum 

We now wish to investigate whether a better fit can be obtained by letting the IC/CMB 
spectrum cut off before the low-energy end of the bandpass under consideration. It is ex- 
pected that if 7min5io ^ 40 (5io = S/10), the cutoff will be too low in energy to affect the 
0.3-7 keV bandpass; obviously, this limiting value of 7min will be higher for the 0.55-7 keV 
bandpass. 

We proceed to fit the spectra with the custom IC /CMB models develop e d in S ection [3l 



For the 0.3-7 keV ba r idpas s, the best fit for both the iBlumenthal fc Gouldl (Il970l ) and the 



Aharonian fc Atoyan] (1l98ll ) kernel, and for both the sharp cutoff in the electron distribution 
or the flat segment case (cases (1) and (2) in Figure Hj), is such that the best-fitting 7„iin is 
above the limiting value. In fact, plotting the fitting statistic as a function of 7min leads to 
Figure [2l which shows that a cutoff due to 7min appears to be detected about 2a statistical 
confidence for two of the four cases. Note that the differences in C-statistic values are 
equivalent to differences in x^, such that the same Ax^ values may be used to determine the 
confidence regions of fitted parameters. 

The best-fit values of the fitting statistic and the constraints that we are able to place 
on 7min are shown in Table [2l In all cases, the fit is as good as the unbroken power law fit 
or better. The y-axis in Figure [2] is normalized to have a A(C-statistic) of for the best 
fit. We expect the previous fit to the unbroken power law to be recovered by pushing 7jnin 
below the limiting value. In detail, the IC/CMB models to not revert exactly to a power law 
even if 'jmin is very low; furthermore, slight differences in the value of the fitting statistic are 
expected near the minimum, given the minimization algorithm employed by XSPEC. The 
limiting value of the fitting statistic as jram 1 can therefore be slightly different from the 
value obtained for the unbroken power law fit above. 

The power law index for the electron distribution is returned in all four cases as s = 
2.6 ± 0.1, consistent with the measurement of the energy index measured for the power 
law fit (a = (s — l)/2). The index is systematically lower for the case of the fiat electron 
distribution, but the difference to the index for the sharp cutoff in the electron distribution 
is not statistically significant. There is almost no degeneracy between the power law index 
and 7min- 

As expected, restricting the analysis to the 0.55-7 keV bandpass eliminates the minimum 
in the fitting statistic, and only upper limits on 7jnin are obtained. The behavior of the fitting 
statistic as a function of 7niin for the restricted bandpass is shown in Figure [31 Again, the 
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relevant data on the best fits and on the 7min constraints are included in Table [2l 

It is worth noting that, contrary to expectations, the upper limits on 7jnin from the 
restricted bandpass are as tight as from the full 0.3-7 keV bandpass. The reason for this 
might be that the residuals to the fit in the restricted bandpass favor an unbroken power 
law down to the lowest bins and are very sensitive to changes in the shape at the low-energy 
end, while the residuals in the bins between 0.3-0.55 keV are more ambivalent about the 
shape of the model in that range. 



4-2.2. Lower Limit on 7min from Optical Flux Measurement 

Figure H] shows the optical and X-ray measurements of PKS 0637-752 as well as the 
predicted IC/CMB spectra based on 7min5io = 60 and the best-fit value of the power law 
index s = 2.6. As 7min becomes smaller, the low-energy extension of the model fitted to the 
X-ray observations w ill approach and at some point over-predict the observed optical flux 



(ISchwartz et a/.ll2000l ). This gives us an opportunity to place lower limits on 7n 



The sum of the individual knots' optical flux density is 0.574 yuJy, and we estimate a 
conservative error on that measurement of 0.1 /iJy. The observed X-ray flux density at 1 
keV is (21 ± 3) nJy. The ratio of optical to X-ray flux densities thus evaluates to 29.0 ± 9.7. 

The IC/CMB models developed in the previous section were adapted to return the 
ratio of optical to X-ray flux as a function of 7min- The effect of the uncertainty on the 
power law index measurement was included by always calculating the minimum ratio over 
the confldence region of the power law index. The calculation of the ratio in the case of 
the flat electron distribution below 7min requires the caveat that the assumption 7 ^ 1 is 
violated, as 70 ~ 3 for energies corresponding to the optical data point. 

In all cases, the ratio is monotonically increasing with decreasing as long as the 
departure from power law behavior occurs at an energy below the reference X-ray energy, i.e. 
when 7min ^lo ^ 75 as determined by the spectral investigation above. For 7min 610 = 50, the 
ratio evaluates to less than 10 and is safely below the observed ratios quoted above. Lower 
limits on 'y^nin are thus obtained by inverting the optical/X-ray ratio vs. 7min relation and 
reading off 7min at the appropriate statistical upper limits on the ratio. The lower limits 
on 7inin thus obtained will actually underestimate the true lower limits, since the two-sided 
errors on the optical to X-ray flux ratio, and not the one-sided upper limits, are used. 

Table [3] summarizes the lower limits on •y^iri thus obtained. Note that at 99% confldence, 
the lower limit in all cases relaxes to the minimum possible value for 7jnin of 1. 
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DISCUSSION 



Tables [2] and [3] summarize the constraints on 7min that we are able to place. If the full 
bandpass of 0.3-7 keV is used, marginal evidence for a break in the spectrum at around 
Tmin'^io = 40-55 emerges, depending on the model for the electron distribution and the 
kernel for IC scattering. Given the uncertainties with the low-energy calibration of Chandra, 
the bins below 0.55 keV might want to be excluded from analysis, in which case only upper 
limits can be placed on 7min from the X-ray spectrum alone: 7niin ^ 80. 

A spurious break due to 7niin might be detected if the column of Galactic absorption 
used in the fit is an under-estimate. However, riu would need to be increased to about 
1.4 X 10^^ cm~^ (an increase of over 50%) to make the fitting statistic increase monotonically 
with 7min- The behavior of tt-h in th e vicinity of PKS 0637 -752 was investigated by visually 
inspecting the IRAS 100 fim map (IWheelock et a/.lll994l ). No evidence for any Galactic 
molecular clouds was found that would average o ut in the large-beam radio surveys resulting 
in a biased value for the quasar's tt-h reported in lDickey &: LockmanI (jl990l ). The difference 
bet ween their value and the one reported in the Lei den/ Argentine/Bonn Survey of Galactic 
HI (IKalberla et all l2005l ) is only 15% (the former being the higher value); it is therefore 
unlikely that uncertainties in the nn measurement alone are responsible for the minimum in 
the fitting statistic. 

In Figure m it can be seen that, for the same value of 7min, the point at which th e spec- 
trum deviates from power law shape is at a lower energv for the lAharonian fc AtovanI (jl98ll ) 
kernel for inverse Compton scattering compared to the iBlumenthal fc Gouldl (119701) kernel. 
by abo ut a factor of 3. This results in tighter lower limits on 7inin for the lAharonian &: Atoyan 
(jl98ll ) kernel, but correspondingly less-strict upper limits. 



The lAharonian fc Atoyaru (Il98ll ) kernel is the more appropriate one for the situation 
of a jet moving relativistically through a homogeneous and isotropic distribution of seed 
photons, as is the case for the cosmic microwave background. The comparison to the other 
kernel can serve to illustrate the order of magnitude of the uncertainty introduced by the 
choice of the inverse Compton scattering formalism. 

We only investigate two possible shapes of the cutoff at 7min: a sharp cutoff such that 
there are no electrons below 7niin (case 1 in Figure H]), and an electron distribution that has a 
constant density below 7iiiin (case 2 in the same figure). The real shape of the cutoff is likely 
to be more complicated than in either of these models, but these two distributions bracket 
the range of expected real distributions, where the cutoff is probably less sharp than in case 
1, but the electron distribution does cut off more quickly than to a constant distribution 
below 7min- As expected, if the cutoff is sharp, the fitting statistic increases rapidly as 7min 
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increases, while the increase is slower for the case of the flat distribution below 7min, since 
the X-ray spectrum in this case exhibits a correspondingly milder cutoff. 



With 7min 5 



10 ^ 80, and assuming F 



6, the equipart ition magnetic field stren gth for a 



jet composed of hot electrons and cold protons is ~ 10/iG ( iDermer fc AtoyanI 120041) , and the 
radio spectrum is expected to extend unbroken down into the hundreds of kHz (jSchwartz 

mi). 



Georganopoulos et all (120051 ) . based on the work of iDermer fc AtoyanI (120041 ). calculate 
the kinetic power in the PKS 0637-752 jet as a function of 7min ^ot the two extreme cases 
of a hadronic jet (equal numbers of electrons and [cold] protons) and a purely leptonic 
jet, assuming that IC/CMB is the dominant emission mechanism for the X-rays, and that 
r = (5. Our limit of 7mm ^lo ^ 80 is consistent with the ranges in 7min considered in 
their paper, for both the leptonic and hadronic jet, and allows for a jet power of around 
5 X 10^^ erg s~^ for a hadronic jet an d < 10^^ erg s~\ for a leptonic jet. Estimations of 
the jet po wer in earlier work, such as E^yecchio et al. (2000) (L = 3 x 10"^^ erg s~^ for 
-7^ ,. = lOl.lCelotti Ghisellini fc Chiabergel J200lh (L = 8 x 10^^ erg s'^ for 7min = 10-20), 
Uchiyama et al\ (120051 ) (L = 9 x 10"^^ erg s~^ for 7min = 20) are likely to be too high. 



or 

given that 7min can in fact be higher than the assumed val ues without violating t he X-ray 
spectral data. Note also that the limits 10 < 7min ^ 40 in lUchiyama et all (120051 ) are too 
restrictive in light of our results. 

For the adopted jet angle with respect to the line of sight of 6.4°, the deprojected length 
of the X-ray jet is on the order of 900 kpc. With the estimated magnetic field strength of 10 
/iG, the cooling time (taking both IC/CMB and synchrotron losses into accoun t) for electrons 
with 7^ 8, 500 is larger than the travel time along the jet (Equation Dl in IStawarz et al. 
20041 ). It would be possible therefore to produce these electrons in the central core and 
propagate them at the bulk velocity to the sites of X-ray emission. 

The high contrast between the X-ray knot and the inner jet in PKS 0637-752 argues 
against the presence of a significant population of electrons with Lorentz factors in the range 
from «ilOO-1000 in the inner jet region at the current epoch. Whether the morphology is 
due to modulated jet activity or to localized particle acceleration sites remains unresolved. 



6. CONCLUSION 

As far as we are aware, this work represents the first attempt at placing confidence limits 
on the low-energy electron distribution cutoff in a jet based on a fit to its X-ray spectrum. 
Usually, papers investigating the broad-band spectra of jet sources simply quote a value for 
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7min that happens to make the model pass through the measured data points. It is very hard 



to co nduct observations in the radio that are able to constrain 7min (iGopal-Krishna. Biermann fc Wiita 



2004 ). especially for the kpc-scale jets, where the frequency at which the cutoff due to 7min 
would manifest itself is below the synchrotron self-absorption frequency. Thus, the soft 
X-ray spectrum as well as any measured optical fluxes are the best tools to shed light on the 
behavior of the low-energy end of the electron distribution. 

We find that 7min 5io ^ 80, which is significantly higher than the value of 10-20 that has 
previously been assumed in broad-band spectral modeling of this source. The kinetic power 
requirement is therefore lessened, but the questions of jet composition and X-ray emission 
mechanism are not addressed conclusively with this finding. 

The present work did not consider other proposed emission mechanisms for the X-ray 
spectrum, such as direct synchrotron or synchrotron self-Compton. Based on the above con- 
siderations, however, IC/CMB remains a viable model for the X-ray emission in PKS 0637-752. 
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Table 1. Chandra ACIS-S Observations of PKS 0637-752 Used in this Investigation' 



ObsID Exposure ObsID Exposure 





Time (ks) 




Time (ks) 


472 


5.52 


1062 


0.62 


473 


3.54 


1063 


1.22 


474 


4.64 


62549 


6.25 


475 


4.64 


62550 


5.20 


476 


2.39 


62551 


5.32 


1051 


1.02 


62552 


5.07 


1052 


1.02 


62553 


4.88 


1055 


2.01 


62554 


11.22 


1056 


1.73 


62555 


4.91 


1058 


1.51 


62556 


4.84 


1060 


0.62 


Total 


78.20 



'^More information on the individual 
observations i s ava ilable in Table 1 of 
Chartas et all tood ). 



- 19 - 



Table 2. Constraints on '')m\n Obtained from the X-ray Spectrum 



X-ray IC/CMB kcrneP c best-fit 7min <5io 7min'5io 

bandpass or power law model distr.^ C-statistic (Icr) (99%) 



0.3-7 keV 


power law 






886.3 








BluGou 


case 


(1) 


879.4 


49 ± 2 


49 ± 6 






case 


(2) 


881.5 


57 ± 4 


< 67 




AhaAto 


case 


(1) 


881.8 


63 ± 2 


< 67 






case 


(2) 


880.2 


72 ± 3 


72 ± 10 


0.55-7 keV 


power law 






817.2 








BluGou 


case 


(1) 


818.6 


< 50 


< 54 






case 


(2) 


817.2 


< 55 


< 63 




AhaAto 


case 


(1) 


817.3 


< 67 


< 70 






case 


(2) 


817.3 


< 73 


< 79 



^BluGou = lBlumeiithal fc Gouldl |l970[ ). AhaAto = lAharonian fc AtovaiJ |l98lh 
''Cases (1) and (2) refer to the different electron distributions in Figure SI 



Table 3. Lower Limits on 7min Obtained from Optical and X-ray Flux Measurement 



IC/CMB kernel e 7mm ^lo 7min<5io 

distr.'^ (la) (99%) 



BluGou 


case 


(1) 


> 4.2 


> 1 




case 


(2) 


> 6.1 


> 1 


AhaAto 


case 


(1) 


> 4.8 


> 1 




case 


(2) 


> 7.7 


> 1 



AhaAto 





Blumenthal & Gould 


Aharonian & Atovan 


( 


1981) 



mm. 



^Cases (1) and (2) refer to the different electron 
distributions in Figure HI 
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Fig. 1. — Extraction regions superimposed on one of the exposures of PKS 0637-752 (ObsID 
62554). The inner and outer jet regions are the two rectangular regions straddhng the jet; 
the background was extracted from the eUiptical region while excluding the counts from 
within the larger rectangular box. These extraction regions are representative of the regions 
used for all observations; however, the change in focus for some of them necessitated larger 
boxes for the inner and outer jet, which were then chosen to again just straddle the image 
of the jet. 
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40 



60 



80 



Fig. 2. — Confidence limits on 7min for tlie IC/CMB model fits, obtaine d from the 0.3-7 keV 
X-ray s pectrum. The bold lines are for the lAharonian fc AtoyanI (jl98ll ) kernel, the thin lines 
for the iBlumenthal &: Gould! (Il970l ) kernel. In both cases, the solid line is for the electron 
distribution with the sharp cutoff at 7min (case (1) in Figure S]), and the dashed line for 
the constant segment below 7jnin (case (2) in Figure H]). The three dotted lines mark the 
A(C-statistic) = 1.00, 2.71 and 6.63 levels corresponding to the 68.3, 90 and 99% two-sided 
confidence limits on 7mm- 
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y ■ 

• mm 

Fig. 3. — Confidence limits on 7jnm for the IC/CMB model fits, obtained from the X-ray 
spectrum using the restricted 0.55-7 keV bandpass. The same conventions as for Figure [2] 
apply. 
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Fig. 4. — Predicted flux of the IC/CMB models developed in Section [3l The adjustable 
parameters are set to 7min = 60 and s = 2.6. The inset plot shows the selected electron 
density distributions, which are all described as a power law with slope s above 7min- Note 
that the normalization of the y axis is arbitrary. The two cases with a modiflcation of 
the spectrum at 7niin are characterized by either the step-function cutoff at 7min (1) or the 
constant density below 7min (2). Case 3 results when 7min is chosen sufficiently low such 
that t he cutoff in the photon spec trum moves outside its bandpass. The spectra based 
on the lAharonian fc AtovanI (jl98ll ) kernel are shown as bold lines, the ones based on the 
Blumenthal fc Gould! (jl970l ) kernel as thin lines. The spectra for the two models coincide 
for case (3). Both sets of spectra have been normalized to the observed X-ray flux. The 
optical data point is also plotted and is violated by the models without a cutoff in the 
electron distribution. The allowable range of 7min5 such that the cutoff occurs below the 
X-ray data point, and the optical flux from the IC/CMB model does not over-predict the 
optical detection, is approximately 5-80. 



